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ABSTRACT 

We probe the dust extinction through the foreground disk of the overlapping 
galaxy pair UGC 5041 by analyzing B,I, and H band images. The inclined 
foreground disk of this infrared-selected pair is almost opaque in B at a projected 
distance of ~ 8 kpc. From the images, we estimate directly the area-weighted 
distribution of differential near-IR extinction: it is nearly Gaussian with 
(r/ — Th) = 0.6 and a = 0.27. For a homogenous dust distribution and a Milky 
Way extinction curve, this corresponds to a face-on distribution p(r) with a 
mean of < ry >= 0.34 and ay = 0.15. For a clumpy dust model the optical 
depth estimate increases to < ry >= 0.41 and ay = 0.19. Even though the 
galaxy pair is subject to different selection biases and our analysis is subject 
to different systematics, the result is consistent with existing case studies, 
indicating that < ry >~ 0.3 is generic for late-type spirals near their half-light 
radii. 

We outline how to estimate from p(r) by how much background quasars 
are underreresented, where projected within ~ lOkpc of nearby spirals, such as 
damped Ly-a absorbers or gravitational lenses; from our data we derive a factor 
of two deficit for flux-limited, optical surveys. 


Subject headings: galaxies: spiral — galaxies: individual (UGC 5041) - 
galaxies: ISM — ISM: dust extinction 
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1. Introduction 

The disks of spiral galaxies contain enough dust to affect significantly the UV and 
optical radiation from the stars within the galaxy and from any sources behind them (see 
e.g., Davies and Burstein, 1992, for an overview). The effect of dust extinction on the 
emerging spectral energy distribution at A < 2/jrn can be subtle, even if large quantities 
of dust are present (e.g., Phillipps et al. 1992, Witt et al. 1992, Rix 1995), and depends 
strongly on the relative layering of the stars and the dust, and the small-scale clumping of 
the dust. The impact of dust extinction on both the colors of the galaxy itself and on the 
statistical extinction and reddening of background sources depends also on the small scale 
clumping of the dust (e.g., Witt and Gordon, 1996; Berlind et al. 1997, hereafter B97) and 
on its concentration towards the spiral arms (White and Keel, 1992, hereafter WK92; Keel 
and White 1995). 

Even though galaxy disks (and their associated dust) may cover only a few percent of 
the sky to redshifts z > 1, they affect particularly interesting lines-of-sight. For example, 
the true rate of damped Ly-a systems depends sensitively on the degree to which the 
detectability of the QSO (usually in the optical) is affected by potential dust extinction in 
the absorber (Bartelmann and Loeb, 1996). Similarly, the incidence of multiple gravitational 
images and of lensing magnification caused by disk galaxies may be altered by the presence 
of dust that may prevent the cataloging of the source (Rodrigues-Williams et.al. 1994; 
Boyle et.al. 1988). 

Measuring the optical depths of spiral disks is most straightforward with background 
light sources, because the influence of dust scattering is minimized, leaving only extinction. 
However, background point sources are insufficient to study the spatial structure of dust 
extinction. Keel (1983) and White and Keel (1992) pioneered a technique to map the 
extinction distribution in a galaxy, by analyzing the light distribution of a background 
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galaxy shining through it. Using three pairs of galaxies, they found that in these systems 
the face-on B-band extinction near the spiral arms was tb ~ 0.3, and tb ~ 0.075 — 0.2 in 
the interarm-regions. 

As the results of such an extinction analysis depend on the objects chosen, the radial 
range probed, the available spatial resolution and the wavelengths of observation, it is 
important to increase the number of case studies and employed techniques to arrive at a 
coherent and comprehensive picture of dust extinction in spiral galaxies. Here we present 
an analysis of the projected galaxy pair UGC 5041, following the approach of WK92. The 
forground galaxy is a highly inclined Sc at cz = 8083 km/s, with a background galaxy at 
cz ~ 12,000km/s, whose center is projected at a mean distance of 14.4" (or 8/i^o 1 kpc ) from 
UGC 5041’s center. 

We decided to present another case study, to complement the existing small set of 
studies already in the literature (WK92 and B97), because the particular galaxy pair, the 
observations and analysis differ in several relevant aspects: 

• The background galaxy behind UGC 5041 was found only in H(1.6/xm) images (Gary 
Bernstein, priv. conv.)\ it is completely obscured in the B-band. Hence, the galaxy 
pair is subject to different selection biases than optically selected pairs and has the 
potential to yield an optical depth estimate of ( tb) *2 1. 

• Given the relatively small projected separation ( R pro j ~ lAR exp ) and qualitatively 
high optical depths (see also NGC 3314, Keel 1983), our analysis is based solely on 
near-IR observations (/ and H). 

• From the observations we calculate the area-weighted optical depth distibution, p(r), 
near the half-light radius of the foreground disk. This is the most immediately 
relevant quantity for assessing the extinction statistics of background sources. 
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• We apply the derived p(r) to explore the under-representation of QSOs lensed by 
spiral galaxies and of QSOs with low redshift damped Ly-a absorption. 

The paper is organized as follows. In §2 we describe the data reduction and analysis. 
In §3 we describe the estimates of the optical depth distribution and its errors; in §4 we 
discuss the implications of this estimate and in §5 we present the conclusions. 


2. Observations and Reductions 

2.1. Observations and Initial Reduction 

The H-band (1.6pm) images were taken on September 19, 1996 at the 2.3m Bok 
telescope on Kitt Peak using the 256x256 IRcam with a scale of 0.2"/pixel. Bias and dark 
frames were obtained at the beginning of the night. The data were taken by integrating on 
the object for 60 seconds, shifting to the sky (~ 5' off-set) for 60 seconds, then returning 
to the galaxy and repeating the procedure. Each telescope dither was to a slightly different 
part of the sky near the object. The total integration time on the galaxy pair was 3600 sec 
(60 object-sky image pairs). In the subsequent analysis, eleven of the final images were 
discarded, because crucial portions of the image were corrupted by bad pixels. 

Most of the data reduction was carried out with IRAFf]. The bias and dark frames were 
subtracted from all galaxy and sky frames. Subsequently, the two sky flats bracketing an 
object image were averaged, normalized and used to flatheld the interleahng galaxy image. 
The resulting images were then aligned using the foreground galaxy center and co-added. 

Reduced and flux-calibrated I-band(7800A) and the B-band(4400A) CCD images of 

3 IR.AF is distributed by the National Optical Astronomy Observatories, which are 
operated by AURA, Inc., under cooperative agreement with the National Science Foundation. 
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UGC5041 were kindly provided by Gary Bernstein. They were originally taken on February 
27 and 29, 1992, respectively, at the 2.3m Bok telescope on Kitt Peak. Since the optical 
images had different pixel scales (0.3"/pixel), the H-band image need to be de-magnihed 
and aligned. All photon-noise and flux-calibration errors were propagated throughout 
these steps, but in the end found to be negligible compared to the systematic modeling 
uncertainties discussed below. 

Figure 1 shows the resulting images in B, I and H. It is apparent that in the B-band 
the extinction is high, for the most part obscuring the background galaxy. Only one bright 
region is shining through (in the upper left part of the foreground galaxy), explaining why 
this galaxy pair was never found in optical surveys[]. Qualitatively, the I-band image hints 
at Tj ~ 1, as the dust extinction is obvious, but not overwhelming. The rightmost panel of 
Figure 1 shows that dust extinction effects in the H-band are subtle. However, we do not 
neglect the H-band extinction. 


2.2. Estimates of tj — th 

With the three images aligned, scaled, and flux-calibrated we can proceed to estimate 
the optical depths in the foreground galaxy. The recorded images are composed of light 
from the foreground galaxy and of attenuated light from the background galaxy. The 
equation for radiative transfer in this case is simply (see also WK92) 

if s (x) = I f r (x) + ^A ack (z) e~ Tx{x) ( 1 ) 

where A denotes the wavelength of the observations, and x the image position. Since the 
extincted source is well behind the extincting dust we neglect scattering into the line of 


Tt was found serendipidously by G. Bernstein in a low resolution H-band image. 
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site. From Eq.(l) one can see that as a first step we must subtract the foreground galaxy 
light, J^ ore (a;), at each wavelength in the region where the two galaxies overlap. Since 
the edge-on background galaxy is an almost linear feature, an estimate of the foreground 
galaxy light was made by interpolating the adjacent foreground light across it. This fit was 
then subtracted from the overlapping region, and the associated errors were an important 
source of uncertainty in the final optical depth determination. Indeed, the Hll regions 
and the spiral arms in the foreground galaxy rendered the B-band image too irregular 
to get an accurate estimate for the fit. Therefore, the B-band data was not used in the 
subsequent analysis. The H-band image on the other hand, had little error in subtracting 
out the foreground light, because the foreground galaxy contributed little to the light in 
the overlapping region to begin with. To get an independent error estimate for Jj) bs — J^ ore 
we analyzed the opposite side of the foreground galaxy, where there are no over-lapping 
galaxies, and carried out the same interpolation. The residuals provide us with an accuracy 
estimate of ~ 6% in 1 for the foreground subtraction. 

Next the corrected background light, I% bs (x) — I { ™ e (x), should be divided by an 
estimate of its unextincted value, I ba,ck (x), to solve for T\(x). Although direct estimates 
of J back (o;) are unknown, because we cannot neglect the foreground extinction even in the 
H-band. Fortunately, the intrinsic 1-H color of highly inclined disk galaxies is quite uniform. 
For example, Bernstein et al (1994) find I-H=1.88±0.08 magnitudes for their sample of 21 
spirals, with 2 outlier galaxies that are bluer by 0.4 magnitudes than the mean. Hence, 
we have a very good guess of the background galaxy’s color in absence of the foreground 
galaxy. By dividing the observed I-H flux ratio by the un-extincted, or intrinsic, value, we 
can estimate Tj — t h directly via 

( yback yintr \ 

T^k / jTntr ) • ( 2 ) 

To obtain a well defined map of 77 — r# we created a mask that isolated the image regions 



of significant background galaxy flux in H. The pixel values in the masked ratio-image can 
then be converted into a histogram of 77 — rwhich is shown in Figure 2. Our data span 
a range of ~ 3.5kpc in the plane of the foreground galaxy, but no systematic radial 
variation of the extinction was found. 

Uncertainties in the optical depths were calculated by standard propagation of errors 
(Bevington 1992), accounting for the foreground subtraction and the intrinsic color ratio 
variance. Intrinsic color ratio errors contributed 0.1 to 77 — th, and the error contribution 
from the subtraction depended on J^ bs — /{ ore . The final median error in 77 — r# was found 
to be ±0.13. 


3. Analysis 
3.1. Estimating p(ry) 

Our data provide a direct estimate of the differential extinction probability distribution 
p(Tj — th). For broader applications it would be useful to correct this value to face-on and 
to convert it to estimates of 7y. The radial range( 6.4 to 10.3 kpc) covered by our estimate 
of p (jv) is quite limited, but of particular interest. First, because damped Ly-a absorbers 
have sizes ~ 10 kpc (Briggs et. al. 1989). Second, because the Einstein ring radius for 
distant galaxies is ~ 4 kpc(e.g. Schneider et. al. 1992). As a first step and benchmark, we 
estimate 7y by assuming an extinction curve, r = r(A), similar to the Milky Way’s, (e.g. 
Mathis, 1990). Under this assumption the face-on value of ry is given by 

T v face -° n = a(r I -T H ) cos i, (3) 

where a ~ 2.15 (Mathis 1990). The inclination angle, i, of the foreground galaxy was found 
by fitting isophotes. Assuming an edge-on axis ratio of 0.19 for UGC 5041’s stellar disk we 
find an inclination angle of 75° ± 2°; the inclination error turned out to be negligible in the 
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subsequent analysis. 

A pixcl-to-pixel histogram of the face-on optical extinction, ry, derived under these 
assumptions is shown in Figure 3. The distribution of the face-on ry is approximately 
Gaussian with (ry) ~ 0.34, and a variance of 0.15. 


3.2. Patchy Extinction Models 

In converting Tj — Th to ry we assumed so far that the extinction curve in the 
foreground galaxy was similar to the Milky Way’s. However, the IR extinction was derived 
over an extended region of the foreground galaxy, as each 0.3" pixel subtends 164 pc 
on a side at a distance of lld/i^Mpc. Casual inspection of nearby galaxy images and 
quantitative studies (see e.g., Witt and Gordon, 1996; WK92, B97) make it clear that the 
dust is quite clumpy on much smaller scales. In this section we investigate how much dust 
inhomogeneities affect our estimates of p(ry) by using a simple model to derive an effective 
coefficient, a e ff, for clumpy extinction. 

Following B97 we consider an area A in the foreground galaxy, covered by one pixel, 
within which a fraction / is covered by material with an X times greater extinction curve 
than the rest of the area. The incoming intensity Jo, emerges as J obs , averaged over the area 
A: 

1 rf A 1 r A 

jobs = / J oe -Xr da + / J oe -r da (4) 

A Jo A Jf a 

Carrying out the integration and denoting 

tteff =-« (5) 

Tv 


where, 


r«e,v = - ln(/ obs //o) = - ln(/e~ Xr + (1 - f)e~ T ) 


(6) 
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gives a e g. In Eq. 6 r is the mean V extinction , equal to 0.34, derived in §3.1. If / is 
set to zero in Eq. 6 then the measured, or effective, extinction will be equal to the mean. 
Qualitatively the situation differs most for / ~ 0.5, where we obtain 

0 496 

« eff = W 2.15 = 3.14, (7) 

for a ratio of optical depths, X ~ 2. Very similar estimates for Ty/{ji — r H ) have been 
derived empirically by B97 (their Figure 3). The histogram for ry = a e //(r/ — t h ), the 
clumpy extinction case, is shown by the dotted lines in Figure 2. The mean extinction in 
V has increased from 0.34 to 0.41, with a commensurate increase of the dispersion. Note, 
that the mean has changed by more than the measurement errors. Therefore, the assumed 
extinction curve is the dominant of uncertainty in our analysis. 


4. Discussion 

4.1. Comparison With Other Authors 

Optical depths through galaxy disks have been derived several ways by different authors 
(Keel 1983, Keel and White 1995, James and Puxley 1995, B97). Each analysis had to rely 
on different sets of assumption, such as intrinsic galaxy symmetry (WK92, B97), case B 
recombination ratios for Balmer lines (James and Puxley 1995), and intrinsic near-IR colors 
of galaxies (this paper). It is interesting to explore to which extent the results are impacted 
by the different techniques. 

James and Puxley (1995) found quite high extinction values from analyzing H/3/Ha 
ratios, seen though the center of NGC3314. These values could be due to the small impact 
parameters of the probed lines of sight, or because the extinction intrinsic to the HII regions 
has been under-corrected. 
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WK92 and B97 have estimated extinction by broad-band photometric methods, 
WK92 using CCD data, B97 using data that extend to 2.2/un. Referring all extinction 
values to face-on and the B-band, we find a mean of ( tb) ~ 0.55 for UGC 5041 at 
R ~ lAR exp ~ 0.4i?25 using our clumpy model, which has an effective extinction curve 
similar to B97. Keel and White (1995) list tb = 0.1 — 0.4 at somewhat larger radii, and B97 
find Tb = 0.6 in an inter-arm region at R ~ 0.7/225- B97 also find tb = 1.3 in the dust-lanes 
near the arms, for which detailed case studies of nearby galaxies (e.g., Elmcgreen 1980 and 
Rix and Rieke 1993) have already shown that t# > 1. 

Even though we have analyzed a purely infrared-selected pair of galaxies, and have 
based our analysis on near-IR data, it appears that our results are quite consistent with the 
optical analyses, implying that the basic result of semi-transparent disks is not simply a 
selection effect. 


4.2. Implications For Cosmologically Distant Objects 

As mentioned above, dust in foreground galaxies will dim and redden distant 
background objects, such as QSOs, that appear in close projection. On the other hand, the 
mass concentration associated with the intervening galaxies may brighten the background 
source through gravitational lensing (Hutchings 1990; Bartelmann et. al. 1996). Hence the 
observed incidence rate of close alignments between background sources and foreground 
galaxies depends on the competition between the two effects. 

Here we show how to use p{jv ) to estimate the first effect, focusing on optically selected 
quasars that lie within ~ (5fol0) kpc of galaxies, such as QSOs lensed by spiral galaxies, or 
QSOs with low-redshift damped Ly-a absorbers. The radial range( 6.4 to 10.3) kpc covered 
by our estimate of p(ry) is quite limited, but of particular interest. First, because damped 
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Ly -a absorbers have sizes ~ 10 kpc (Briggs et. al. 1989). We take from Peterson (1997, 
p. 174) the integrated number-magnitude relation for quasars brighter than Mf 

N(< Af t ) = N 0 10 a8Mi = N 0 e L84Mi (8) 


For objects whose light passes through a foreground galaxy, the magnitude limit of any 
survey is effectively M[ = Mi — 1.086r, where r is the optical depth at the wavelength 
of the survey. Given a probability distribution p(r) we would expect such sources to be 
under-represented by a factor of: 


iVobs(< MA 


N{< M,) lV 0 e L84M * Jo 

As Figure 3 shows p{r) is well approximated by 

1 

p(r) = 


N oe ( M i-i-osr) dT ' 


(9) 


(- t ~(t)) 


\/2tT<7 t 


which leads to 


iVobs(< Ml) = -^^(l-^) 
N(< M t ) 


( 10 ) 


( 11 ) 


For the derived values of (r) = 0.34 and oy = 0.15, and (r) = 0.41, o> = 0.19, this results 


in = °- 53 and = 0.47, respectively. 


Hence the chance of finding lensed of damped Ly-o: absorbed quasars (with m B < 19) 
in such close projection (< 10 kpc) to spiral disks is suppressed by about a factor of two 
due to dust extinction, in good agreement with the independent estimate based on the 
Milky Way gas-to-dust ratio by Perna, Loeb and Bartelmann (1997). 


5. Conclusion 

We have presented and analyzed H, 1 and B band images of the projected galaxy pair 
UGC 5041, in order to explore the distribution of dust extinction that light suffers while 
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crossing the highly inclined foreground galaxy disk. This galaxy pair is completely obscured 
in the optical, and hence could have lead to (ly) > 1; this system also provides optical 
depth constraints at quite small impact parameters. 

After removing the light contribution from the foreground galaxy, we compared, pixel 
by pixel, the observed I-H color of the (extincted) background galaxy to its expected 
intrinsic 1-H color, which can be estimated to ±0.08 mag. This comparison resulted in 
a direct estimate of the distribution of 77 — ry, with (77 — th)/ cos 72° = 0.6. For a 
homogeneous dust distribution this translates into a face-on V-band optical depth of 
(tv) = 0.34, for clumpy dust this corresponds to (tv) = 0.41, with dispersion of 0.15 and 
0.20 in Ty. 

Despite the different selection criteria (which would have permitted to arrive at much 
higher optical depths) and an analysis based solely on near-IR data, these optical depths 
are consistent with the values found by WK92 and B97, and corroborate that generically 
^face— on ^ j n near by galaxy disks. 

Using the derived distribution of p(ry), we estimated by how much quasars in close 
projection (5 — 10 kpc) are under-represented in optically selected, flux-limited samples 
(rrihmit(B) < 19), and found a deficit by a factor of two due to dust extinction. 

As in many other analyses, the limited spatial resolution, which leads to averaging over 
areas of differing extinction is an important limitation. We are currently analyzing HST 
images of UGC5041 to repeat an analogous analysis at vastly improved spatial resolution. 

We would like to thank Gary Bernstein for providing the CCD frames UGC5041. HWR 
was supported by an Alfred P. Sloan Fellowship. 



-14- 

REFERENCES 

Berlind, A., Quillen, A., Pogge, R. and Sellgren, K., 1997, AJ, 114, 107 (B97). 

Bernstein, G. et al. , 1994, AJ, 107, 1962. 

Bartelmann, M., Loeb, A., 1996, ApJ, 457, 529. 

Bevington, P., Robinson, D.,1992 in Data Reduction And Error Analysis for the Physical 
Sciences, ed. S.J.Tubb,J.M.Morriss(McGraw-Hill,Inc.)p.43. 

Boyle, B. J., Shanks, T., 1988, MNRAS, 231, 897. 

Briggs, F., Wolfe, A., Liszt, H., Davis, M. and Turner, K., 1989 ApJ 341 650. 

Davies, J., Burstein, D. in The Opacity of Spiral Disks, Kluwer Academic Publishers, 1995. 
Elmegreen, D., 1980, ApJS, 43, 37. 

Hutchings, J.B., 1990, PASP, 102, 431. 

James, P., Puxley, P.1995, in The Opacity of Spiral Disks, ed. J.I.Davies, D.Burstcin(Kluwer 
Academic Publishers)p.l81. 

Keel, W.C., 1983, AJ, 86, 1579. 

Keel, W.C., White, R.E., 1995, in The Opacity of Spiral Disks , eds. J.I.Davies,D.Burstein, 
Kluwer, p.167. 

Mathis, J.S., 1990, ARAA, 28, 37. 

Perna, R., Loeb, A. and Bartelmann, M., 1997, ApJ, 488, 550. 

Peterson, B.M., 1997, in An Introduction to Active Galactic Nuclei, (Cambridge University 
Press, p.159. 

Phillipps, S., Evans, R., Davies, J. Disney, M., 1992, 253, 496. 

Rix, H.-W. and Rieke, M., 1993, ApJ, 418, 123. 



- 15 - 

Rix, H.-W., 1994, in The Opacity of Spiral Disks, eds. J.I.Davies, D. Burstein, Kluwer, 
p.215. 

Rodrigues-Williams, L.L., Hogan, C., 1994, AJ, 107(2), 451. 

Schneider, P., Ehlers, J. and Falco, E., 1992, in Gravitational Sensing , Springer Verlag. 
White, R., and Keel, W.C., 1992, Nature, 359, 129 (WK92). 

Witt, A.N., Thronson, H.A., Capuano, J.M., 1992, ApJ, 393, 611. 

Witt, A.N. and Gordon, K., 1996, ApJ, 436, 681. 


This manuscript was prepared with the AAS IAT^X macros v4.0. 



16 





Fig. 1.— The B, I, and H band images of the galaxy pair UGC5041, from top to bottom. 
The center of the two galaxies are separated by 14.4". The images cover 39" x 64", and 
35" x 66", and 38" x 53" in B, I, and H, respectively. 







Fig. 2.— The pixel-by-pixcl histogram of the differential optical depth 77 — r#. The error 
(±0.124) is discussed in Section 2.2. The distribution can be approximated by a Gaussian 


with < tt — tw >= 0.6 and oy = 0.27. 




Number 


18 





Fig. 3.— Histogram of the face-on optical depth in V, assuming a homogeneous dust 
distribution (solid line) and a clumpy dust model (dashed line). The errors are 0.07 and 
0.08, respectively. 























